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Using classical differential geometry, the problem of elastic curves and surfaces in the 
presence of long-range interactions <I>, is posed. Starting from a variational principle, the 
balance of elastic forces and the corresponding projections • V3>, are found. In the case 
of elastic surfaces, a force coupling the mean curvature with the external potential, K^, 
appears; it is also present in the shape equation along the normal principal in the case 
of curves. The potential ^ contributes to the effective tension of curves and surfaces and 
also to the orbital torque. The confinement of a curve on a surface is also addressed, in 
such a case, the potential contributes to the normal force through the terms —K^ — n-^^. 
In general, the equation of motion becomes integro-differential that must be numerically 
solved. 
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1. Introduction 

The elastic energy of a polymer chain, can be modeled as a geometric functional, in- 
variant under reparametrizations, H = J ds /(k, t), k, and t being the Frenet-Serret 
curvatures and s the arclength of the curve accounting the polymer. Working within 
the geometric framework of Frenet-Serret, an infinitesimal deformation of the curve 
can be decomposed in its perpendicular and tangential parts; the invariance under 
reparametrizations has the effect that just orthogonal projections play a role in the 
equilibrium shape equations while a tangential deformation yields only a boundary 
term. Is possible to integrate the Euler-Lagrange equations for the curvature in 
the case of plane curves If the energy only depends on the curvature, the model 
is integrable, and then the torsion is a function only of curvaturtPl. Furthermore, 
Noether conserved charges of the model, associated with invariances of motions in 
space, can be obtained using this geometric formalism. 

In certain processes in cell biology and in certain technological applications, 
the behavior of polymers under the influence of external fields is of great interest. 
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The problem of polymers chains constrained to confined geometries was already 



on the polymer. For example, in the case of DNA packaging, not only the entropic 
and elastic forces act on the polymer, but also the electrostatic force. Indeed, as the 
DNA is immersed in a solution with several other ions, the effective electrostatic 
potential can be modified as the Debye-Hiickel potential SI cx r~^e~^^ ^ that 
results from a linearized mean-field procedure. In general all polyelectrolytes, which 
are charged polymer chains, that are in solution have similar behavior. When the 
polymer chains are negatively charged, the charges produce their mutual repulsion 
and the polymer can not be wrapped, so the chain is stretched They also have 
certain behavior when the chain is close to a charged surface, this behavior has 
been used to model some biological interactions such as in nucleosome structures 
of DNaI^. The properties of the polyelectrolytes are well known, see e.g. refs. [THSl 
however, the geometric properties under external long-range interactions are not 
yet studied in detail. In the present work we consider besides the elastic forces, the 
long-range interactions on curves and surfaces. 

In section ([2]), starting from a variational principle, we obtain the Euler-Lagrange 
equations for elastic curves in presence of long range forces. As a result of transla- 
tion and rotational invariance of the energy, the corresponding Noether charges are 
found. The example of elastic curves dependent on the curvature k, and the par- 
ticular case of planar curves are shown, for these cases a first integral of the shape 
equation can be found. In section ([3]), the analogous analysis for surfaces with long 
range forces is presented. The shape equation from a variational principle and the 
corresponding forces and torque are found. Explicit equations in the case of fluid 
membranes are presented. In section we use the alternative method introducing 
auxiliary variables to regain the Euler-Lagrange equations of section The equa- 
tions of elastic curves in presence of long range forces, constrained to surfaces are 
presented. Moreover, the normal force that constraint the curve is found in terms 
of the geometric information. We conclude with some final remarks in section ([S]). 

2. Curves 

2.1. Shape equations 

In the presence of long range interactions, an elastic curve can be modeled by the 
following energy functional 




[ However, there was not taken into account the role of a potential acting 




(1) 



The first term Ha = / dsf{K,T), constitutes the elastic energy of a curve in R^, 
X = X(s), parametrized by arclength s^. The curvature k and torsion t of the 
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curve are defined through the Frenet-Serret equations: 

t = Kill, 

hi = -Kt + rn2, (2) 
h2 = -mi, 

here the dot indicates tlie derivative respect to the arclength s. The unit tangent 
vector t = X, together the unit normal vectors n.;, define a local basis {t, rii, n2}, 
adapted to the curve^i^. Among others, some relevant cases are the energy /(k) = 
associated with the bending of a curve, the Euler elasticc^^, the Sadowski f unc tional 
/(k) = (k^ + T^Y / hi^ related with bending of elastic developable surfaces^^, etc. 

The second term in the model ([T|) is Tib = 5// dsds'h{r) with r = |X(s) — 
X(s')|, represents the energy of the curve in the presence of long range forces, 
electrostatic interaction /i(7') = l/r as well as the screened Coulomb interaction 
(Yukawa) h{r) = e^°'^'/r and the Lennard- Jones potential h{r) = (7)^^ — (^)^ 
some relevant possibilities to explore. The invariance under euclidean motions as 
well as reparametrizations of the model can be exploited to obtain several relevant 
properties. 

In order to obtain the shape equations we have to consider an infinitesimal de- 
formation of the energy 51-1, as a result of deformations in the embedding functions, 
X — >■ X + i5X. For instance, given that the arclength is defined by ds^ = dX • dX, 
its deformation can be written as Sds = d((5X) • t. Thus, and following reference [131 
we can see that an arbitrary infinitesimal deformation can be written as 

SHa = J ds 5X ■ni + J ds Q, (3) 

where £i (i = 1,2), are the Euler- Lagrange operators and Q the corresponding 
Noether charge. 

Let us now consider, the deformation of the second term in equation ([T]). Since it 
is symmetric under the change of s and s', we have that SHb = J S{ds) J ds'h{r) + 
/ ds (5X • V / ds'h{r), where V denotes the 3D gradient operator. After integration 
by parts and introducing the potential $ := J ds' h{r) we can write 

SHb = J d{SX-t<P) - J dsKSX-ni<^>- J dsSX-fl> + J dsSX-\/<P. (4) 

Notice that the last term of equation ^ can be decomposed in its normal and 
tangential projections as 5X ■ V$ = SX • t <i> + {6X ■ nj)ni • W^. The first term 
cancels the corresponding tangential projection in equation Q. As expected by 
invariance of Ji under reparametrizations, tangential deformations do not play a 
role in the shape equations, these can be obtained from SH = 6Ha + SUb = 0: 

£;i := 5i - K$ + ni • V$ = 0, 

E2 :- £2 + n2 • V$ = 0. (5) 

The normal projections of the external force V$, have contribution to these equa- 
tions as it should be. It is interesting to note that the potential $, contributes 
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directly to the shape equation along the principal normal direction ni. This might 
be interpreted as a minimal coupling between the curve and the external field. 

2.2. Translational invariance 

By adding equations ([3|) and (|4]), the total deformation of the energy H can be cast 
as 

Sn^ j ds Ei (5X ■ni+ J dsQ (6) 

where now, Q = Q + (5X • t <!>. is the corresponding Noether charge. According 
to equation ([5]), under a translation (5X ~ a. the energy is deformed as SH — 
a. - J ds (^Ei rii — Gj , where G = F — t <&, being the relation with the elastic forces 
through Q = a ■ F. In equilibrium G is a conserved vector field along the curve. 
Expressing in the local basis, G = G||t + GiUi, the balance equation EiUi = G 
yields 

G|| - kGi = 0, 
Gi+kG|| -rG2 =^1, (7) 

G2+TG1 ^E2, 

the first of them, or Bianchi identity, implies that the arclength of the curve is 
preserved under the total force G. In terms of the elastic and external forces we 
have 

F\\ - kFi = t • V$, 
Fi + - TF2 = £1 + ni • V$, (8) 
F2 + tFi =£2 + n2- V$. 

The first equation tell us that arclength is not preserved under the elastic force F, 
but it is balanced with the tangential component of the external force V$. In the 
same way, this force contributes to the equilibrium in the normal directions as we 
can see from the last two equations in ([5]). 

2.3. Rotational invariance 

The invariance under rotations can be explored in a similar way as in ref. 1131 
Considering an infinitesimal rotation (5X = x X in equation ([6|), we obtain 

sn = n- I ds (e, X X n, - m) , (9) 



here, the relation Q ~ — Jl • M, was used. The equation i?i X x = M. follows from 
rotational invariance. M is conserved along curves in equilibrium and it is identified 
with the torque respect to the origin. The long range potential contributes to the 
orbital torque as 

M = Xx (F-t$) + T, (10) 
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where the local torque T satisfies that T = G x t. Moreover, T = F x t since 
the contribution of the external force is in the direction of the tangent vector t. 
Therefore, decomposing T in the local basis, it satisfies 

T|| - kTi = 0, 

ti-TT2 + nT^\^F2, (11) 

Ta+rTi = -Fi, 
the external long range potential $, does not contribute. 

2.4. Elastic curves 

Let us focus on the case where the elastic density energy is only function of the cur- 
vature / = /(k). The two corresponding Euler-Lagrange operators and the Noether 
charge are given by 

Ei=h + {k.^ - T^)U - ^^{f + $) + ni • V$, (12) 

E2 = 2tA + f + n2 • V$, (13) 

Q = (/ + $)*|| + fJi - f^^i - 2tU^2, (14) 

where = df /dn and the projections of the deformation, = ^||t + ^'in^, were 
used. The conserved force G, is given by 

G= (/.'«-/-$)t + /Kni+r/„n2, (15) 

thus, the external potential <I> contributes to the tension in the curve. There is no 
contribution of the external force to the local torque T and then T = —fK,'n-2 • 
Furthermore, we can write a first integral of the problem, in terms of the conserved 
quantities G and J = G • T: 

G^ = /.' + (/«^-/-$)' + ^- (16) 

The presence of the long range force, into the effective potential V{k, $) — (/kK — 
/ — $)^, does not allow their interpretation as a central potential V(k) as in the 
elastic case. 



2.5. The Euler elastica 

In the case of plane curves with bending energy /(k) = k^/2, the Euler-Lagrange 
equation reduces to: 

K + K - $^ = -n- V4>, (17) 

this was found in a dynamical context in ref. 1141 The corresponding first order 
integral equation can be written as 

«^+(y-$)' = G^ (18) 



January 11, 2013 2:7 WSPC/INSTRUCTION FILE Curves-and-Surfaces- 
R 



6 J. A. Santiago, G. Chacdn-Acosta and O. Gonzdlez-Gaxiola 

and for the particular case of electrostatic interaction it is given by 

in terms of the embedding function X(s), it is an integro-differential equation. 

3. Surfaces 

3.1. Shape equation 

The geometric description of the surface is constructed considering that it is em- 
bedded in R"^. We parameterized it by local coordinates a ~ 1,2 , through 
X = X(^°). Now there are two tangent vectors Ga = SqX to the surface, and the 
unit normal vector field n is defined by Gq • n = and n • n = 1. The induced metric 
on the surface is defined by the symmetric tensor gab = • g^, whose inverse is 
denoted g"-^. We denote Vq the covariant derivative compatible with the induced 
metric. 

The model of the surface that we consider is the sum of two terms, a natural 
generalization of equation ([T]) for curves 

^ = y" dAf{K) + ^ J J dAdA'h{r), (20) 

where dA = ^ d^^ is the infinitesimal area element on the surface, g being the 
determinant of the induced metric. 

The first term in equation ([20)) . J^i, is the elastic energy of the surface; f{K) is a 
scalar function constructed with the geometry of the surface. K = g°-^Kab the mean 
curvature in terms of the second fundamental form Kab = — VaG;, • n. The second 
term in equation (j2Gp . Jif2, involves the external potential h(r), similar to the case 
of curves, section ([2]). The functional Jif, is invariant under reparameterizations and 
euclidean motions of the surface. 

Given the invariance under reparameterizations of the ener gy, o nly normal de- 
formations (5X • n play a role in the Euler- Lagrange equations Deforming the 
first term of the energy we can obtain 

(5^1= J dA£SX-n + J dAS/aQ"", (21) 

where £ is the Euler-Lagrange derivative of the elastic model and the corre- 
sponding Noether charges. The remaining deformation of the energy J^f in eq. ([20|l 
involves the external potential. We obtain 

S-M ^ j '^^ (-^^'* + " • 5^-^ + j dA\Ia {g^^SX ■ Gf,*) , (22) 

where we have defined $ = / dA' h{r). The shape equation can be found by setting 
5Jfi + 5,yif2 = 0. From equations ([2T|l and ([22| we have that the Euler-Lagrange 
derivative is 



(23) 
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in equilibrium, the elastic forces are balanced with the external long range force, 
£ = — /1<1> — n • V$. In the case of minimal surfaces where K ~ 0, there is no 
contribution of the second term to the shape equation. 



3.2. Translational and rotational invariance 

We can write the total deformation of the energy in the form 

S,y/^^ J dA{£ + /v$ + n • V$) 5X • n + J dAVa (Q" + g'^'^SX ■ efa$) . (24) 

By doing an infinitesimal translation (5X = a we obtain that 

SJ^f ^a- J dA[(£ + A'$ + n • V$) n - (f ° - e° $)] , (25) 

where the relation between the Noether charge and the elastic forces, = —a - f", 
was used. Thus, on a surface in equilibrium, h° = f ° — e° $ is a conserved vector 
field, as a consequence of invariance under translations. Using the decomposition in 
the local basis, = /"''eh + /"n, we can write 

Va/" - KabP'' - £ + n • V$, (26) 

comparing with eq.® we can see that they match for the case of curves. 

Let us see the consequences of the invariance of the energy under rotations. Let 
(5X = n X X be an infinitesimal rotation in the deformation of the energy equation 
(|24)) . We have 

(5^ = n • y (i^[^X X n - VaM"] , (27) 

the invariance under rotations of the energy implies that (fX x n = VaM°. Thus, 
in equilibrium M" is a conserved vector field on the surface and it is identified with 
the torque. We can decompose this vector in an orbital torque plus a local one, in 
the form M° = X x (P — e"$) + s" where the local torque satisfies VaS" = x Gq. 



3.3. Lipid membranes 

In this casJi^ the energy is quadratic in the mean curvature f{K) = The 
Euler-Lagrange derivative and the conserved Noether charge of the elastic energy 
are modified as follows 

= -\/^K + K (^2Kg - ^K^^ + + n • V$ 
^ K {2K''^ ~ Kg''^) ef, - 2V° A'n + $ e" (28) 
where 2Kg — TZ, is the curvature scalar of the surface. 
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4. Constraining curves on surfaces. Auxiliary variables. 

It is possible to impose the relations ([2]) directly in the energy ([1]) by considering 
them as constraints of the variational problem, an d th en fixing them through the 
introduction of indeterminate Lagrange multipliers In this way, the curvature 
K and X can be varied independently. In order to confine the curve into a surface 
X = Y(^") we have to add in equation ([T)) this constraint in the form 

Hc^H + J dsX-[:S.{s)-Yie{s))]. (29) 

As in ref. |2|we take the variation with respect to X obtaining the following relation 

SxHc = J ds (f + X + V$) • SX, (30) 

where V$ denotes the euclidean gradient, and the multiplier F is given by 

F ^ ^ nf^)t + (-fr + h) ni + I rf^ - Kfr - ^ I — ]] n2, (31) 




At being the Lagrange multiplier that enforces the curve to be parameterized by its 
arclength. From equation we see that along a curve in equilibrium, F + V$ = 
—A. We can also see that A is orthogonal to the surface and therefore we can 
write, A = —An. As a consequence the tangential projection vanishes identically, 
i.e. ^F + • Bq = 0, in particular, projection along the tangent t determines the 
multiplier At in equation ([?T|l to be: 

At = ./ + $ - 2k/k - rfr + a, (32) 

where cr is a constant that tell us that the length of the curve is fixed. Thus, the 
external potential contributes to the effective tension. Projections of F + V<I>, 
along the normal vectors n^ give the Euler-Lagrange equations. We have 




(F + V$) • ni = + (k^ - t')U + 2t— ( i± ) + r ( ^ I + 2Kr/, 

- K(/ + $ + cr) + ni • V$, (33) 

(F + V$)-n2 =2-^(T/J-f/,- 
as 

' ■"'+n2-V$, (34) 




K 

these equations reproduce the Ei derivatives of equations ([S]) in the case of curves 
without restriction. Then we can write (F + V$) = £^11^. If we identify the normal 
to the surface to be the principal normal of the curve, n = ni, then we can write 



F + V$ = EnU + Ell, 



(35) 
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where the unit vector field 1 = n2 is defined by the conditions 1 • t = and 1 • n = 0. 
From equation ([55]) . the Euler Lagrange equation, _Ei = and the normal force, —A, 
can be recognized as 



Ei = {F + V$) • 1, (36) 
-A = (F + V<I>) • n. (37) 

That is, the long range force contributes to the normal geometric forces adding the 
terms k$ — n • V$. 

With this formulation one can also consider the long range interaction between 
points on the surface with points on the curve, this can be achieved by adding the 
term 

J J dsdA'h{R), (38) 

where now R = |X(s) — Y (^°(s')) |. This term considers the contributions to the 
interaction of the surface on the curve, for every point on the curve. Notice that the 
back-reaction that the curve has on the surface is not considered. The force satisfies 

F + V (^<I> + = -A(s), (39) 

where $ = / dA' h{R) is the contribution of the potential on the surface. Bearing 
in mind that now we can vary with respect to the coordinates on the surface 
that leads us to 

S^n = - J ds A(s) • eaS^" - J ds Vl> • eaS^", (40) 

such that (A(s) + V$) ■ Ba = 0. This means that the only force that appears on the 
surface is V<i>. 

From the variation of H with respect to the remaining variables one can obtain 
the Lagrange multipliers: 

A«=H«, Ar=nr. (41) 

Using the properties of stationarity and Frenet-Serret equations, the following re- 
lation for the components of the multiplier F, which corresponds to the force, can 
be derived: 

F = (-A - kH^) t+(-nr+ n.) ni + (tH^ -nUr-^ (k-^Ht) ] na, (42) 



K 



ds 



where only the multiplier A is missing. To find A one can take the derivative of (|^^ , 
and compare its components with those obtained in p9l) , which are clearly different 
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form zero. With this, one can find the following relations: 

F-t = -(K-H^)--rH^-KH«-A = -$, (43) 
F • ni - {ti^ + T^) Tin + 2r (k^^H^^ + f (k^^H^) + KrHr - nX 

= - (^V$ + A + V$) -ni, (44) 

F • ns = {rn^y- (k-^H,)"- {K-HrY+T {j4Lr+H^) 

= - fv$ + A + Vl>) -nz. (45) 



From (|43p one can integrate to obtain A, noticing that this expression is a total 
derivative on s: 

X = H- 2kH^ -rUr + '^ + cr, (46) 

where a again fixes the length of the curve. With the multiplier F is determined 
as: 

F = {kH^ +THr-H-^-(T)t + {^Ur + Hk) Hi 

tU,- KUr- (K-^nr))n2. (47) 

Amusingly, both the force that keeps the curve on the surface and the interaction 
with the membrane, do not appear in this expression; only the potential due to the 
interaction of the curve itself appears in the tangential component. 

5. Concluding Remarks 

The interest on studying charged polymers and surfaces comes because in nature, 
there are some systems, such as the polyelectrolytes El, that can be modeled as a 
charged chains which require a description involving long-range electrostatic forces. 

In this work we have studied curves and surfaces under the influence of long- 
range and elastic forces, both from the classical differential geometry perspective. 
Besides the usual elastic terms, we introduce the long-range interactions as an inte- 
gral of the corresponding potential, which takes into account the interaction between 
all points of the curve ([1]). The corresponding shape equations were obtained form 
a variational principle, considering small deformations. 

In the case of curves the resulting Euler-Lagrange equations can be written as 
the sum of the usual elastic term and normal projections of the external force. 
Amusingly, in one normal directions there are a contribution of the bare potential, 
which further is coupled with the curvature. We can interpret this term in ([5|) as a 
minimal coupling. Given the translational and rotational invariance of the model, 
we were able to write the equations for the force and torque, noticing that the 
long-range potential $ contributes directly to the orbital torque, unlike the force 
and the local torque, where it only contributes through its gradient. In the case 
where the elastic force is quadratic on the curvature, the well known Euler elastica 
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model, the corresponding equation for k becomes an integro-differential equation 
that wiU be solved numerically elsewhere. For the charged surfaces we perform the 
same analysis straightforwardly, and found a similar minimal coupling, now with the 
mean curvature, and the projection of the long-range force in the normal direction; 
both in the Euler-Lagrange derivative. 

We also consider the case of a curve constrained to lie on a surface, this analysis 
was done by introducing Lagrange multipliers for the constraints of the system, 
recovering the previous results, and obtaining additional information for the normal 
force on the curve. Moreover, with this formalism, we could write the corresponding 
interaction in the case where both the surface and the curve were charged. 

To integrate all the resulting equations, numerical methods must be used, al- 
though it was not the main task of this work, we are in the process of implementing 
this in a subsequent study. 
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